THE CONTRIBUTIONS OF INTERACTIVE BINARY STARS 
TO DOUBLE MAIN SEQUENCE TURN-OFFS AND DUAL 
RED CLUMP OF INTERMEDIATE-AGE STAR CLUSTERS 



Wuming Yang^'^, Xiangcun Meng^, Shaolan Bi^, Zhijia Tian^, Tanda Li^, Kang Liu^ 

ABSTRACT 

Double or extended main-sequence turn-offs (DMSTOs) and dual red clump 
(RC) were observed in intermediate-age clusters, such as in NGC 1846 and 419. 
the DMSTOs are interpreted as that the cluster has two distinct stellar popula- 
tions with differences in age of about 200-300 Myr but with the same metallicity. 
The dual RC is interpreted as a result of a prolonged star formation. Using a 
stellar population-synthesis method, we calculated the evolutions of binary-star 
stellar population (BSP). We found that binary interactions and merging can 
reproduce the dual RC in the color-magnitude diagrams of an intermediate-age 
cluster, whereas in actuality only a single population exists. Moreover, the binary 
interactions can lead to an extended MSTO rather than DMSTOs. However, the 
rest of main sequence, subgiant branch and first giant branch are hardly spread 
by the binary interactions. Part of the observed dual RC and extended MSTO 
may be the results of binary interactions and merger. 

Subject headings: galaxies: star clusters: general — globular clusters: general — 
binaries: general 

1. Introduction 

It has always been accepted that a star cluster (SC) is composed of stars belonging to 
a single, simple stellar population with a uniform age and chemical composition. In recent 
years, however, due to the increase in spatial resolution and photometric accuracy, it has 
been discovered that some SCs have unusual structures in their observed color-magnitude 
diagrams (CMDs). For example, NGC 2173 in the Large Magellanic Cloud (LMC) has 
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an un usually large spread in color about its main-sequence turn-o ff (MSTO) (IBertelli et al. 



20031 ): Omega Centauri exhibits a main sequence (M S) bifurcation (iPiotto et al.ll2005l ): NGC 



pie stellar populations ( iMilone et al. 



2808 possesses a triple MS split (IPiotto et al.ll2007D : and other clusters also contain muti- 



2009 



20101). More interesting discoveries are that the 



intermediate-age massive clusters in the LMC, such as NGC 18 46, 1806, and 1783, pos- 
sess double main-sequen ce turn-offs (DMST Os) on thei r CMDs JMackev fc Brobv Nielsen! 
200?! : iMackev et allEoOsi iMilone et allboooh . Moreover. Icoudfrooii et all ( l2009h found that 



the sp read of the MSTO is f airly c ontinuous rather than strongly bimodal. iMackey et al. 



( 120081 ) and iGoudfrooij et al.l (120091 ) argued that these clusters do not possess a significant 
line-of-si gth depth or inter n al dis persion in [Fe/H], or suffer from significant differential ex- 



tinction. 



MucciareUi et al.l fl2008f ) showed that spreads in [Fe/H] in SCs like NGC 1783 are 
very small. The apparent homogeneity in [Fe/H] of all stars in NGC 1846 indicates that cap- 
ture of pre-existing field stars during cluster formation seems to hard to explain the DMSTOs 
( iMackey et al.ll2008t IGoudfrooij et al.ll2009l ). The DMSTOs are interpreted as that the clus- 
ters have two distinct stell ar populations with differences i n age of ~ 200-300 M yr but with 
similar metal abu ndance ( iMackey fc Broby NielsenI 120071 : IMackey et al.l 120081 ) . Moreover, 
Glatt et al.N2008l) found multiple MSTOs (MMSTOs) in NGC 419 in the Small Magellanic 
Cloud, and iGirardi et al.l (l2009[ l found two distinct red clumps (RC) in this cluster. These 
observational results are challenging the traditional picture. 

In order to understand how a clust er can contain multiple popula tions with differ- 
ent ages, many scenarios were proposed. IMackey fc Brobv NielsenI (120071) pu t forward the 



Goudfrooij et al. 



(120091) pointed out 



scenario of merger of two (or mor e) SCs. However , ^ ^ _^ ^ 

that this scenario seems unlikely. Bekki fc Mackey ( 20091 ) simulated the merger of a SC 
with a giant molecular cloud (CMC). They found that this merger can produce a second- 
generation of stars that are required to explain the DMSTOs. However, in order to obtain 
the differences in age of ~ 200 - 300 Myr, this scenario require rather strongly constrained 
ranges of G MC parameters such as their spatial distribution, m ass function, and che mical 



composition (iGoudfrooii et al.ll2009l : iBastian fc de MinkI 120091 ). iD'Ercole et al.l (120081 ) and 



Goudfrooij et al. ( 20091 ) put forward the scenario of formation of a second generation of stars 



from the ejecta of first generation asymptotic giant branch stars. In this scena rio, the young 



second generation is usually l ess than the first on e ( IBastian &: de MinkI |2009| ). which is not 



consistent with observations ( Milone et al.l 2009 ). Rec ently, using two scaling relations to 



mimic the effects of rotation, IBastian &: de MinkI ( 120091 ) considered the effect of rotation on 



stellar evolutions. They found that stellar rotation in stars with masses between 1.2 and 
1.7 Mq can mimic the effect of a double population when rotation rates reach 20-50 per 



Rubele et al. 


(2010 


) and 


Girardi et al. 


(2011) 
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the presence of MMSTOs. They proposed a prolonge d star formation (PSF) to explain 
the MMSTOs and the dual RC. In iRubele et al.l (120101 ) models, in order to reproduce the 
dual RC of N GC 419, a PSF of a bout 700 Myr is required. The dual RC was also found 
in NGC 1751 jRubele et allbojA In order to explain this dual RC, a PSF of about 460 



Myr is re quired (IRube. 



mated bvlMilone et al. 



eetal 



201ll ). This age spread is twice longer than ~ 200 Myr esti- 



._ j2009h using the method of isochrone fitting. Furthermore, recently, 

Keller et al. ( 201ll ) argued that the isolation of extended MSTO clusters to intermediate 
ages could be the consequence of observational selection effects. 

In addition, binary s tars of LMC clusters are unresolved even with Hubble Space Tele- 
scope (IMilone et al.ll2009l ). The binary system should appear as a single point-source object. 
Thus the unresolved binary systems may have c ontributions to the spr e ad of MSTO. How- 
ever, several investigators (IMackey et al.l l2008l : iGoudfrooij et al.l 120091 : iMilone et al.l 120091 ) 
had found that unresolved binaries cannot alone reproduce the peculiar MSTO structures. 
Nevertheless, after a star in a binary system accretes mass from its companion, it will be- 
come you nger apparently. Hence, interactive binaries may have contributions to the spread 
of MSTO. IGoudfrooij et al.l ( l2009l ) found that the younger populations in clusters are more 
centrally concentrated than the older populations. Thus the interactive binaries may be 
important in interpreting the observed CMDs. In this p aper, using the Hurley rapid single 
and binary evolution codes (IHurley et al. I l2000l . |2002|) and a stellar population-synthesis 
method, we investigated the contributions of interactive binaries to the DMSTOs and dual 
RC. 

The paper is organized as follow. We show our stellar population-synthesis method in 
section 2. We present the results in section 3 and discuss and summarize them in section 4. 



2. Stellar population synthesis 

To investigate DMSTOs and dual RC of SCs, we calculated a binary-star stellar pop- 
ulation (BSP). Stellar samples are generated by Monte Carlo simulations. The basic as- 
su mptions f o r the simulations are as follows, (i) The lognormal initial mass function (IMF) 



of IChabrierl ( 1200 ll ) is adopted, (ii) We generate the mass of the primary. Mi, according to 
the IMF. The ratio (g) of the mass of the secondary to that of the primary is assumed to 
be a uniform distribution within 0-1 for simplicity. The mass of the secondary star is then 
determined by qMi. (iii) We assume that all stars are members of binary systems and that 
the distribution of separations (a) is constant in log a for wide binaries and falls off smoothly 
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an{a) 



a. 



(1) 



at close separation: 

t/ao)™ a < ao] 
sep, ao < a < fli, 

where ttsep ~ 0.070, oq = IORq, ai = 5.75 x lO^i?0 = 0.13pc and m ^ 1.2. This distribution 
imphes that the numbers of wide binary system per logarithmic interval are equal, and that 
approxi mately 50% of th e stellar systems are binary systems with orbital periods less than 
100 yr ( Han et aD Il995l ). However, binaries account for typicall y ^30-40% of all 
stars in the cluster s of LMC, such as in NGC 1818 and 1806 (.Elson et al.lll998l : 



Milone et al.l |2009| ) . (iv) The eccentricity (e) of each binary system is assumed to be a 



uniform distribution within 0-1. With these assumptions, we calculated the evolutions of 
5x10^ binaries with Mi between 0.8 and 5.0 Mq. 

The metal abundance Z of evolutionary models was converted firstly into [Fe/H]. Then 
the theoretical properties ([Fe/H], Te//, iogg, logL) have been transformed into colors and 



magnitudes using the color transformation tables of iLejeune et al. I (119981 ). The binaries 
with a < 10^ R.Q were treated as unresolved ones, while others were treated as resolved ones 
when we computed their colors and magnitudes. 



3. Calculation results 



Figure [T] shows the CMDs of the simulated BSP with Z = 0.008 and age = 1.8 Gyr. 
The interactive binaries are shown in red, while non-interactive binaries are shown in green. 
The population of the brighter MSTO (bMSTO) is mainly from the unresolved binaries with 
g > 0.7 and no interactions. The isochrones of single-star stellar populations (SSP) with the 
same Z but different ages are overplotted on the CMD in the right panel of Fig. [H The stars 
that experienced binary interactions clearly deviate from the isochrone of non-interactive 
binaries in the region around MSTO, which leads to a large spread in color. The interactive 
binaries clearly produce an extended MSTO rather than bimodal MSTO. For main sequence 
binary stars with my < 22.5 mag, about 10% of the binary stars lie between the 
1.8 and 1.5 Gyr isochrones in our models, not including the binaries on the 
1.8 Gyr isochrone. However, interactive binaries hardly affect the rest of MS, subgiant 
branch and first giant branch (FGB). Moreover, interactive binaries result in the appearance 
of a secondary RC (SRC) at about 0.5 mag below the main RC. The SRC is slightly bluer 
than the ridgeline of FGB stars. The contribution of non-interactive binaries to the SRC is 
negligible. The right panel of Fig. [1] shows that the ridge of the SRC is almost coincident 
with the isochrone of RC stars of the SSP with age = 1.2 Gyr. About 12% of core-helium 
burning (CHeB) binary stars are located between the 1.8 and 1.5 Gyr isochrones. However, 
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about 30% of CHeB binary stars are located between the main RC and 1.2 Gyr isoclirone, 
not including the binary stars on the main RC but including those on 1.2 Gyr isochrone. 
The mass of the SRC stars is mainly located between 1.85 and 2.1 Mq, while that of the 
main RC stars is less than 1.7 Mq. 



The dual R C was discovered in NGC 752, 7789 and 419 flGirardi et all l2000l 120091 ) 



and NGC 1751 flRubele et all l201lh. Cluster NGC 419 has an age of about 1.35 Gyr 
(ICirardi et al.l 120091 ). We computed a cluster with age = 1.3 Gyr and Z = 0.004. Its CMDs 
are shown in Fig. |21 Just as above results, the interactive binaries lead to a large spread 
in color in the region around MSTO, but they do not reproduce DMSTOs. However, the 
interactive binaries lead to the presence of a distinct SRC. The ridge of the SRC is almost 
coincident with the isochrone of RC stars of the SSP with age = 1.1 Gyr. This implies that 
the SRC produced by binary interactions has an apparent age of about 1.1 Gyr in our mod- 
els, whereas actually it has an age of 1.3 Gyr. For this dual RC, about 15% of binary 
stars are located between the 1.3 and 1.1 Gyr isochrones, including the binaries on the 
1.1 Gyr isochrone but not including those on 1.3 Gyr isochrone. The magnitude extension 
between main RC and SRC is about 0.35 mag. The mass of the SRC stars is mainly located 
between 1.85 and 2.05 Mq, while that of the main RC stars is mainly located between 1.78 
and 1.81 Mq. For Z = 0.004, in stars more massive than about 1.88 Mq (this value is 
about 2.0 Mq for Z = 0.02), helium-burning temperatures are reached at the center before 
electrons become degenerate there. For these stars, when they are on zero-age horizontal 
branch (ZAHB), their magnitudes decrease with decreasing mass. However, for the stars 
with M < 1.88 Mq, electrons in the hydrogen-exhausted core are highly degenerate before 
helium ignition occurs. When they are on ZAHB, their magnitudes increase rapidly with 
decreasing mass except for stars with M < 1.5 Mq. The magnitudes of ZAHB stars with 
M < 1.5 Mq increase very slowly with decreasing mass. For Z = 0.004, the evolutionary 
tracks of RC stars with mass between about 1.85 and 2.0 Mq are very near each other in 
CMDs. Thus the interactive binaries with mass between about 1.85 and 2.0 Mq can gather 
to form a SRC. Our calculations show that the SRC caused by binary interactions appears 
mainly in clusters with 1.2 Gyr < age < 3.0 Gyr. 



Discussion and Conclusions 



Milone et al.l (120091 ) showed that the fraction of younger (brighter) population is about 



75% in the case of NGC 1846. In our models, the 'younger' population are mainly from the 
merged binary systems and interactive binaries with mass transfer. The merging and mass 
transfer in binary systems can be affected by the distributions of separation (a), eccentricity 
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e) and mas s-ratio (q) of system s. We computed the BSP by using the distr ibutions n flog a 



= constant flHurley et al. II2002I ). n{e) = 2e and n{q) = 2q and the IMF of ISalpeterl ( 119551 ) . 
However, the fraction of the interactive binaries are not obviously enhanced. We also com- 
puted the BSP using the distribution r;,(loga) = constant and Gaussian distributions for e 
and q. The mean value and standard deviation of the Gaussian distributions is 0.5 and 0.13 
for e, 0.6 and 0.1 for q, respectively. The values of the standard deviation are chosen in order 
to make the values of e and q are located between and 1. In this simulation, about 20% 
of binary stars lie between the 1.8 and 1.5 Gyr isochrones for the cluster with 
Z = 0.008 and age = 1.8 Gyr; while for the RC of the cluster with Z = 0.004 
and age = 1.3 Gyr, about 30% of binary stars lie between the 1.3 and 1.1 Gyr 
isochrones. 

The mass of turn-off stars is about 1.46 Mq for the cluster with Z = 0.008 and age = 
1.8 Gyr in our models. The evolutions of a wide binary system without mass transfer are 
similar to those of single stars with the same mass. However, the evolution of stars whose 
mass increased via mass accretion or merging is always slower than that of single stars with 
the same mass. For example, a mass- accreted star with M = 1.55 Mq, which was reached 
by accretion, just arrives at around MSTO at the age of 1.8 Gyr. However, a single star 
with an initial mass of 1.55 M© has left the MSTO at the same age. Thus the mass-accreted 
stars clearly deviate from the ridge of the isochrone in the region around MSTO and have 
an apparent younger age. The more the accreted mass, the more the deviation. Thus the 
spread of MSTO of BSP is continuous rather than bimodal. However, for the mass-accreted 
stars with M < 1.25 Mq, an increase in mass leads to a shift of their position in CMDs 
almost along MS. Therefore, the most of MS are hardly spread by binary interactions. 

The SRC is almost made of the mass-accreted or merged stars whose mass is slightly 
more massive than the critical mass that is just enough to avoid electron-degeneracy occur- 
ring in stellar H-exhausted cores in our models. The stars with the critical mass usually 
evolve into RC at the age of about 1.1-1.2 Gyr (For the stars with Z > 0.04, this age can 
increase to about 1.3 Gyr). Hence, the SRC caused by interactive binaries should have an 
apparent age of about 1.1-1.2 Gyr, whereas actually it has an age as old as that of main 



RC. The ap parent age of about 1.1-1.2 Gyr is consistent with the estimates of iRubele et al. 



(j2011 . 2010) for NGC 1751 and 419. The main RC stars passed through electron degeneracy. 



The magnitude of these stars increases with decreasing mass. The older the age of clusters, 
the lower the mass of main RC stars. Thus, the magnitude extension between main RC and 
SRC increases with age until the mass of the main RC is less than about 1.5 Mq. Because 
the magnitude of ZAHB stars with M < 1.5 Mq increases very slowly with decreasing mass, 
there is an upper limit of magnitude extension of about 0.8 mag between main RC and SRC 
in our models. Moreover, with increasing age, the SRC becomes slightly bluer than FGB. 
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For the clusters with age < 1.2 Gyr, the mass of all RC stars is almost larger than the critical 
mass. Therefore, the SRC caused by interactive binaries should not appear in these clusters. 
In addition, in old clusters, the initial mass of RC stars is low. If the mass is not enough 
to be increased to above the critical mass by accretion or merging, the SRC should not be 
produced by binary interactions in the old cl usters. Moreo v er, if the fraction of binaries is 



too low, the SRC also could not be produced. iGirardi et al.l (120091 ) argued that the MS+red 



clump binaries cannot mimic the dual RC in NGC 419. In our models, the SRC stars are 
from the merged binaries and binary systems with mass transfer. 

In this paper, we showed that binary interactions such as mass transfer and binary 
merging can produce an extended MSTO and dual RC in CMDs of intermediate-age clusters, 
whereas in actuality only a single population exists. Despite these, the rest of MS, subgiant 
branch and FGB are not clearly spread by the binary interactions. Interactive binaries can 
lead to an extension of MSTO rather than DMSTOs. For a cluster with Z = 0.008 and age 
= 1.8 Gyr, its isochrone can be spread down to ~ 1.5 Gyr by interactive binaries. However, 
only about 10% of binary stars lie between the 1.8 and 1.5 Gyr isochrones in our 
models. Moreover, the SRC that is caused by interactive binaries should have an apparent 
age of about 1.1-1.2 Gyr. Although binary interactions cannot lead to the bimodal MSTO, 
the SRC in NGC 419 and extended MSTO of NGC 1846 may be partly from the mass- 
accreted or merged binary stars. 
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Fig. 1. — The CMDs of simulated BSP. Green shows the binaries almost without mass 
transfer, while red indicates interactive binaries. A distance modulus of 18.45 is adopted. 
The isochrones of SSP with the same Z but with different ages are overplotted on the CMD 
in the right panel (blue: 1.8, cyan: 1.5 and magenta: 1.2 Gyr). 
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Fig. 2. — Same as Fig. [H A distance modulus of 18.84 is adopted. The isochrones of SSP 
with the same Z but different ages are overplotted on the CMD in the right panel (blue: 1.3 
and cyan: 1.1 Gyr). 



